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Abstract —Microbial communities play a significant role in 
bioremediation, plant growth, human and animal digestion, 
global elemental cycles including the carbon-cycle, and water 
treatment. They are also posed to be the engines of renewable 
energy via microbial fuel cells which can reverse the process 
of electrosynthesis. Microbial communication regulates many 
virulence mechanisms used by bacteria. Thus, it is of fundamental 
importance to understand interactions in microbial communities 
and to develop predictive tools that help control them, in order 
to aid the design of systems exploiting bacterial capabilities. This 
position paper explores how abstractions from communications, 
networking and information theory can play a role in understand¬ 
ing and modeling bacterial interactions. In particular, two forms 
of Interactions in bacterial systems will be examined: electron 
transfer and quorum sensing. While the diffusion of chemical 
signals has been heavily studied, electron transfer occurring in 
living cells and its role in cell-cell interaction is less understood. 
Recent experimental observations open up new frontiers in the 
design of microbial systems based on electron transfer, which may 
coexist with the more well-known Interaction strategies based 
on molecular diffusion. In quorum sensing, the concentration of 
certain signature chemical compounds emitted by the bacteria is 
used to estimate the bacterial population size, so as to activate 
collective behaviors. In this position paper, queuing models 
for electron transfer are summarized and adapted to provide 
new models for quorum sensing. These models are stochastic, 
and thus capture the inherent randomness exhibited by cell 
colonies in nature. It is shown that queuing models allow the 
characterization of the state of a single cell as a function of 
interactions with other cells and the environment, thus enabling 
the construction of an information theoretic framework, while 
being amenable to complexity reduction using methods based on 
statistical physics and wireless network design. 

Index Terms —Quorum sensing, electron transfer, bacterial 
interactions, stochastic modeling, queuing models 

I. Introduction 

Bacteria constitute some of the earliest life forms on Earth, 
existing anywhere from four to three billion years ago ||T]. 
Bacteria are unicellular organisms that lack a nucleus and 
rarely harbor membrane-bound organelles. While bacteria have 
certain elements in common (such as being unicellular), their 
numbers and variety are vast, with an aggregate biomass 
larger than that of all animals and plants combined Q. 
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Bacteria exist all over the planet, in and on living creatures, 
underground, and underwater. They can survive in Antarctic 
lakes and in hot springs, showing tremendous robustness and 
tenacity. It is theorized that bacteria enabled changes in the 
Earth’s environment leading to atmospheric oxygen as well as 
being the precursors to more complex organisms. Microbial 
communities play a significant role in bioremediation, plant 
growth, human and animal digestion, global elemental cycles 
including the carbon-cycle, and water treatment |[^. Despite 
their simplicity as unicellular organisms, the operations and 
interactions of bacterial communities are not fully understood. 
Elucidating how bacteria interact with each other and with 
their environment is of fundamental importance in order to 
fully exploit their potential. To this end, realistic tools for the 
modeling and prediction of these systems are needed 

In this position paper, we explore how queuing models can 
play a role in understanding and modeling bacterial inter¬ 
actions. Indeed, molecules, electrons and cells are countable 
units. It is thus natural to represent them as "quanta", and their 
amount as the state of a "queue" in which these quanta are 
collected and from which they are dropped. Queue evolution 
then depends on the interaction of these elements with each 
other and with the surrounding environment. The advantage of 
a queuing model compared, for instance, to continuous models 
based on ordinary differential equations, is their amenability 
to capture randomness and interactions in small bacterial 
systems, such as the quorum sensing system studied in Q, 
where the discrete nature of these interactions is predominant. 
In particular, we shall examine two forms of interactions; 

a) Electron Transfer & While chemical and 

molecular diffusion are widely investigated electron 

transfer is emerging as an exciting strategy by which bacteria 
exchange nutrients and potentially information. Each cell relies 
on a continuous flow of electrons from an electron donor to 
an electron acceptor through the cell’s electron transport chain 
to produce energy in the form of the molecule adenosine 
triphosphate (ATP), and to sustain its vital operations and 
functions. This strategy, known as oxidative phosphorylation, 
is employed by all respiratory microorganisms. While the 
importance of biological electron transport is well-known 
for individual cells, the past decade has also brought about 
remarkable discoveries of multicellular microbial communities 
that transfer electrons between cells and across much larger 
length scales than previously thought Q, spanning from 
molecular assemblies known as bacterial nanowires, to entire 
macroscopic architectures, such as biofilms and multicellular 
bacterial cables 0, ®. Electron transfer has been observed 
in nature 0 and in colonies cultured in the laboratory m- 
This multicellular interaction is typically initiated by a lack of 
either electron donor or acceptor, which in turn triggers gene 
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expression. The overall goal in these extreme conditions is for 
the colony to survive despite deprivation, by relaying electrons 
from the donor to the acceptor to support the electron transport 
chain in each cell. The survival of the whole system relies on 
this division of labor, with the intermediate cells operating as 
relays of electrons to coordinate this collective response to the 
spatial separation of electron donor and acceptor. 

These observations open up new frontiers in the modeling 
and design of microbial systems based on electron transfer, 
which may coexist with the more well-known interaction 
strategies based on molecular diffusion. A logical question to 
ask is whether classical approaches to cascaded channels E). 
GD are suitable for the study of bacterial cables, which are 
indeed multi-hopped networks (T3). A straightforward cut-set 
bound analysis 0 shows that the capacity of such a system 
is achieved by a decode-and-forward strategy. However, in 
microbial systems, there is a stronger interaction between the 
bacteria, due to the coupling of the electron signal with the 
energetic state of the cells via the electron transport chain. 
Thus, alternative strategies for analysis and optimization must 
be undertaken^ which take into account the specificity of 
the interactions between cells in a bacterial cable. In turn, it 
is important to understand how to design capacity-achieving 
signaling schemes over bacterial cables, which exploit these 
specific interactions, as a preliminary step towards optimizing 
fuel cell design (see Sec. |I-B| i. In Sec. [I^ we summarize a 
stochastic queuing model of electron transfer for a single cell, 
which links electron transfer to the energetic state of the cell 
{e.g., ATP concentration). This model abstracts the detailed 
working of a single cell by using interconnected queues to 
model signal interactions. In Sec. Ill we show how the pro¬ 
posed single-cell model can be extended to larger communities 
{e.g., cables, biofilms), by allowing electron transfer between 
neighboring cells. In Sec. |III-A| based on such a model, we 
perform a capacity analysis for a bacterial cable, and discuss 
the design of signaling schemes. 


b) Quorum Sensing 0-0-- As previously noted, 
biological systems are known to communicate by diffusing 
chemical signals in the surrounding medium ||7|, Q. In quo¬ 
rum sensing, the concentration of certain signature chemical 
compounds emitted by the bacteria is used to estimate the 
bacterial population size, so as to regulate collective gene 
expression. However, this simplistic explanation does not fully 
capture the complexity and variety of quorum sensing related 
processes. In this paper, we provide a model for quorum 
sensing, informed by our work on electron transfer, that will 
enable more sophisticated study and analyses. 

While quorum sensing is based on the emission, diffusion 
and detection of these chemical signatures across the environ¬ 
ment in which the cell colony lives, it differs from recent work 
endeavoring to model and design transceiver algorithms based 
on molecular diffusion channels (see, e.g., 0-0) and to 
analyze their capacity (see, e.g., p^-p6)). The focus of these 
works is the design of models and algorithms for engineered 
systems for future nanomachines to act as transmitters and 


'Relay systems employing decode and forward for molecular diffusion 
based communication systems are considered in 0 


receivers. Herein, instead, we attempt to understand natural 
systems and their optimization via the control of environmental 
conditions. Rather than looking at the characteristics of the 
diffusion channel, we consider the system (the cell colony and 
the surrounding environment) as a whole, and attempt to model 
the dynamic interactions between cells and of the cells with the 
environment. We abstract the diffusion channel by considering 
a queueing model to represent it in terms of concentration of 
signals. As will be seen in our model depicted in Figs.j^an d0 
quorum sensing does not fit well into traditional multi-terminal 
frameworks such as the multiple access channel p7]-|30|, the 
broadcast channel 1(31), 0 or two-way communication 0. 
In fact, recent work examining capacity questions relating to 
quorum sensing have considered the binding of molecules and 


multicellular processes with molecular diffusion |25|, 

The models adopted therein do not take the dynamics of the 
binding processes into consideration as we do here. 


A. Prior Art on Mathematical Modeling of Bacterial Popula¬ 
tions 

The modeling of natural populations such as swarms (motile 
bacteria, swimming fish, flying birds, or migratory herds 
of animals) has been consistently studied over the years. 
Microscopic models model individuals as point particles via 
dynamical equations These models have precision, 

but become much more complex as the interactions are finely 
modeled and do not scale well as the population size increases. 
To combat such issues, the computational approach of simulat¬ 
ing a large number of agents with prescribed interaction rules 
has also been considered 0, 0. These studies duplicate 
experimental results well, but do not lead to design methodolo¬ 
gies or optimization strategies as the underlying operations are 
not considered. In contrast, macroscopic models model biolog¬ 
ical swarm dynamics via advection-diffusion-reaction partial 
differential equations 0, 0-0. The partial differential 
equations are used to describe the evolution of the probability 
distribution of group sizes or fraction of individuals with 
specific characteristics {e.g., motion orientation |[42)). These 
approaches also suffer from computational complexity, often 
relying on numerical solvers, and do not have an eye towards 
optimization and control. 

From an information theoretic perspective, there have been 
ongoing efforts to examine how entropy and mutual informa- 
tior0 can provide insights into experimental data |[44)-||4^; 
however these works do not endeavor to model the underlying 
processes, but rather compute empirical information theoretic 
measures based on experimental data. Another key distinction 
is a lack of an explicit modeling of the dynamics of the 
constituent systems. 

Our modeling approach differs from these mathematical 
models on the following aspects: 1) We propose an accurate 
queueing approach to characterize the state of a single cell as 
a function of interactions with other cells and the environment. 
Thus, not only this approach can be used to describe the 
metabolic state of a cell, but it also enables the construction of 
an information theoretic framework; 2) our proposed models 

^See |43| for a review of the application of information theory to biology. 
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BACTERIAL CABLE 


Fig. 1. Stochastic model of electron transfer within a bacterial cable. 

are inherently stochastic (rather than deterministic), and thus 
capture the inherent randomness exhibited by cell colonies in 
nature; 3) finally, they are amenable to complexity reduction, 
e.g, using methods based on statistical physics gz) or wireless 
network design | [4^ . We show an example of application of 
these principles to the capacity analysis of bacterial cables in 
Sec. HIkAl 

B. Applications 

The proposed queueing models serve as powerful predictive 
tools, which will aid the design of systems exploiting bacterial 
capabilities. We highlight two relevant applications. 

1) Microbial Fuel Cells: Renewable energy technologies 
based on bioelectrochemical systems are now attracting tens 
of millions of dollars in government and industry funding Q. 
Microbial fuel cells and the essentially reverse process of 
microbial electrosynthesis are well known bioelectrochemi¬ 
cal technologies; both use microorganisms to either generate 
electricity or biofuels in a sustainable way p^ . Microbial 
fuel cells are constructed with microbes that oxidize diverse 
organic fuels, including waste products and raw sewage, 
while routing the resulting electrons to large-area electrodes 
where they are harvested as electricity pO). Certain bacterial 
strains (e.g., Shewanella oneidensis lig, |51| ) are of particular 
interest for microbial fuel cells as they can attach to electrodes 
and transfer electrons without mediators which are often toxic, 
especially in the high concentrations needed to overcome their 
diffusion. 

The challenge to realizing microbial fuel cells is under¬ 
standing the cooperative and anti-cooperative behaviors of 
collections of bacteria. Ideally, as the number of bacteria 
increases, the amount of current generated should increase as 
well. However, this is not always observed in many experi¬ 
mental scenarios fS^ . We posit that the cell-cell interactions 
necessitated by a multicellular (i.e. biofilm) lifestyle, including 
electron transfer across and through the biofilms themselves, 
can limit device performance. In other words, the biofilm’s 
energy output cannot be casually gleaned from what we know 
about the mechanisms relevant to individual or “disconnected” 
microbes. Optimization of biofilms, and multicellular perfor¬ 
mance in general, requires sophisticated models and predictive 
tools that account for the electronic and nutrient fluxes within 
whole communities. In this paper, we propose queuing models 
to serve this purpose. 


2) Suppression of Infections: Quorum sensing processes 
have been of significant interest since first observed in Vibrio 
harveyi |Tg. The discovery of quorum sensing in Pseu¬ 
domonas aeruginosa, a human pathogen associated with cystic 
fibrosis and burn wound infections, demonstrates the medi¬ 
cal implications of quorum sensing regulation | |53| , | |54) . In 
Pseudomonas aeruginosa infections, quorum sensing plays an 
important role in activating many critical virulence pathways 


|56|. The genes regulated by quorum sensing enable 


the bacteria to mount a successful attack or evade the body’s 
defenses. In these cases, blocking quorum sensing activation 
may help to prevent or treat microbial-associated diseases gz)- 
As bacteria build resistance to antibiotics, reducing virulence 
through inhibition of quorum sensing may be more effective 
in mitigating disease versus eradicating the bacteria | [58) . 
Naturally occurring compounds and enzymes exist which can 
inhibit quorum sensing gT), However, the complexity 
of signaling interactions within complex communities of mi¬ 
crobes points to the need for a more quantitative understanding 
of quorum sensing regulation. Since many different types of 
bacteria use a similar set of signals, and many signals work in 
conjunction with other autoinducers within the same cell, the 
potential outcomes of quorum sensing manipulation are not 
obvious. The combination of predictive model development 
with quantitative experiments to parametrize and test models 
are essential to develop successful strategies for virulence 
reduction. 

This paper is organized as follows. In Secs. 0 and [m] 
we propose a model of electron transfer for a single cell 
and a bacterial cable, respectively, and present some capacity 
results. In Sec. |IV[ we present a model of quorum sensing 
in a homogeneous bacterial community. In Sec. |V] we present 


some simulation results and experimental data and, in Sec. VI 


we discuss extensions of the proposed quorum sensing model. 
Finally, in Sec. |VII[ we conclude the paper. 


H. Single Cell Model oe Electron Transeer 


In this section, we propose a queuing model of electron 
transfer in a single cell. We will extend this model to inter¬ 
connections of cells in a bacterial cable in Sec. Im] The cell 
is modeled as a system with an input electron flow coming 
either from the electron donor via molecular diffusion, or 
from a neighboring cell via electron transfer, and an output 
flow of electrons leaving either toward the electron acceptor 
via molecular diffusion, or toward the next cell in the cable 
via electron transfer (Fig. [T]i. Inside the cell, the conventional 
pathway of electron flow, enabled by the presence of the 
electron donor and the acceptor, is as follows (see the numbers 
in Fig.[T]): electron donors permeate inside the cell via molec¬ 
ular diffusion (1), resulting in reactions that produce electron- 
containing carriers (e.g., NADH) (2), which are collected 
in the internal electron carrier pool. The electron carriers 
diffusively transfer electrons to the electron transport chain, 
which are then discarded by either a soluble and internalized 
electron acceptor (e.g., molecular Oxygen) or are transferred 
through the periplasm to the outer membrane and deposited 
on an extracellular electron acceptor (3). The electron flow 
through the electron transport chain results in the production 
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TABLE I 

Description of the queuing abstraction for the electron 

TRANSFER MODEL. 


Queue 

Arrival rate 

Service rate 

Explanation 


'^F.XT 

m 

Mryt 

High energy external membrane 


AL) 

'^F.XT 

{L) 

i^FXT 

Low energy external membrane 

■mo Hip 

^OH 


Internal electron carrier pool 

nATp{t) 


IdATP 

ATP pool 


^IextW S 2 


(H) 


EXT — {Oj 1) ■ • • ) ^^Exri^) S 




Q 


(i) 

EXT 

{L) 


= { 0 , 1 ,..., 


Hi 

}, respectively, where Q^^xt 


Qext '■h® respective electron “storage capacities”. The 
internal cell state at time t is thus 

si{t) = (mcff(0,nATp(i),gBirW>9£LxVw) • (1) 

Cell behavior is also influenced by the concentrations (Td( 1) 
and aE{t) of the electron donor and acceptor in the surround¬ 
ing medium, respectively. Therefore, we define the external 
state as SE{t) = {aD{t),(TA{t)). 

Each pool in this model has a corresponding inflow and 
outflow of electrons, each modeled as a Poisson process with 

^While in the following analysis we assume that one queuing "unit” 
corresponds to one electron, this can be generalized to the case where one 
"unit" corresponds to Ne electrons. 

^Note that the terms high and low referred to the energy of electrons are 
used here only in relative terms, i.e., relative to the redox potential at the cell 
surface. In bacterial cables, the redox potential slowly decreases along the 
cable, thus inducing a net electron flow. 


of a proton concentration gradient (proton motive force, @) 
across the inner membrane of the cell (4), which is utilized 
by the inner membrane protein ATP synthase to produce ATP, 
collected in the ATP pool and later used by the cell as an 
energy source (5). 

A. Stochastic cell model 

In this section, we present our proposed cell model. Each 
cell incorporates four pools (see Eig. [^1, each with an associ¬ 
ated state, as a function of time t: 

1) The internal electron carrier pool, which contains the 
electron carrier molecules (e.g., NADH) produced as a result 
of electron donor diffusion throughout the cell membrane 
and chemical processes occurring inside the cell, with state 
mcnit) G Mch = {0, 1 ,..., Men} representing the num¬ 
ber of electron^ bonded to electron carriers in the internal 
electron carrier pool, where Mch is the electro-chemical 
storage capacity, 

2) The ATP pool, containing the ATP molecules produced 
via electron transport chain, with state nATp{t) G J^axp = 
jo, 1, ..., Naxp}, representing the number of ATP molecules 
within the cell, where Naxp is the ATP capacity of the cell; 

3) The external membrane pool, which involves the extracel¬ 
lular respiratory pathway of the cell in the outer membrane. 
This is further divided into two parts; a) a high energj^ 
external membrane, which contains high energy electrons 
coming from previous cells in the cable; and b) a low en¬ 
ergy external membrane, which collects low energy electrons 
from the electron transport chain, before they are transferred 
to a neighboring cell. We denote the number of electrons 
in the high energy and low energy external membranes as 



Fig. 2. Stochastic model for an isolated cell after the transient phase, and 
Markov chain with the corresponding transitions, for the case where Mch = 
^AXP = 4. The transition rates from state (2, 2) are also depicted. 

rate function of the internal and external state {sj{t),SE{t)). 
The intensities of these Poisson processes are labeled as in 
Pig. 0 Por notational simplicity, the input and output flows are 
denoted as A and p,, respectively, which commonly correspond 
to arrival and service rates in the queueing literature. A 
description of the queuing model abstraction is given in 
Table H 

B. Model Validation for an Isolated Cell 

The experimental investigation of a multi-cellular network 
of bacteria is very challenging, due to the technical difficulties 
of placing multiple cells next to each other in a controlled way, 
and of performing measurements of electron transfer, ATP and 
NADH concentrations along the cable 0 . Therefore, we start 
by investigating the properties of single, isolated cells, which 
constitute the building blocks of more complex multicellular 
systems. 

In the case of an isolated cell, multicellular electronjransfer 
does not occur, hence P^ext ~ ^^ext ~ ^ 1 for an 

explanation of the intensities Pext ^ext^ Therefore, 
after a transient phase during which the low and high energy 
external membranes get emptied and filled, respectively, the 
cell’s state and the corresponding Markov chain and state 
transitions are as depicted in Pig.|^ In |j^, we have validated 
this model based on experimental data available in | [62) , using 
the following parametric model for the rates, inspired by 
biological constraints (e.g.. Pick’s law of diffusion J^): 

Acir(s/(f);SE(f)) =p(l- 

POUT{^l{t); SE{t)) = PCH{si(t);SE{t)) 




= C 1- 


nATP jt) 
Naxp 


CtAit), 


( 2 ) 


PATpisi{t)-SE(t)) = / 3 aD(,t), 


where p,C,,fS G M+ are parameters, estimated via curve fitting, 
and Mch — Naxp — 20. 

Initially, cells are starved and the electron donor concen¬ 
tration is 0. At time t = 80s, glucose is added to the cell 
culture. Afterwards, the electron donor concentration profile 
decreases in a staircase fashion from aoit) = 30[mM] 
glucose at time t = 80s, to aoit) = 0[mM] at time 
t — 1300s, and aoit) — 0[mM] for t < 80s, whereas 
the electron acceptor concentration (molecular Oxygen) is 
constant throughout the experiment, and sufficient to sustain 
reduction (aAit) = 1, Vf) ^2\. We have designed a parameter 
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Fig. 3. Prediction of expected ATP level over time (in mM) and experimental 
time-series HD 

estimation algorithm based on maximum likelihood principles, 
outlined in to fit the parameter vector x = [p, ^,/3]. 

Fig. plots the ATP time-series, related to the cell culture, 
and the expected predicted values based on our proposed 
stochastic model. We observe a good fit of the prediction 
curves to the experimental ones. In accordance with |j^, the 
addition of the electron donor to the suspension of starved cells 
triggers an increase in ATP and NADH production as well as 
ATP consumption, as observed experimentally and predicted 
by our model. Further discussion can be found in | |6T| . 

III. Towards a Model of Bacterial Cables 

In multicellular structures, such as bacterial cables, an 
additional pathway of electron flow may co-exist, termed 
intercellular electron transfer, which involves a transfer of 
electrons between neighboring cells, as opposed to molecules 
(electron donor and acceptor) diffusing through the cell mem¬ 
brane. In this case, one or both of the electron donor and 
the acceptor are replaced by neighboring cells in a network 
of inter-connected cells. This cooperative strategy creates a 
multicellular electron transport chain that utilizes intercellular 
electron transfer to distribute electrons throughout an entire 
bacterial network. Electrons originate from the electron donor 
localized on one end of the network and terminate to the 
electron acceptor on the other end. The collective electron 
transport through this network provides energy for all cells 
involved to sustain their vital operations. 

Since typical values for transfer rates between electron 
carriers {e.g. outer-membrane cytochromes) on the cell exterior 
are relatively high | [l3) , we assume that = oo 

when two cells are connected, so that any electron collected 
in the low energy external membrane is instantaneously trans¬ 
ferred to the high energy external membrane of the neighbor¬ 
ing cell in the cable. Therefore, the low energy and high energy 
external membrane pools of two neighboring cells can be 
joined into a common pool for intercellular electron transfer. 
The internal state of cell i at time t is thus 

sW(t) = , (3) 

and its external state is {f)={a'^^ {t), {t)). 

The state of the cable of length N is denoted as 


s(f)=(s(i)(f),..., s(^)(f)), where s^^(f)). 

The internal state Sj\t) is time-varying and stochastic, 
and evolves as a consequence of electro/chemical reactions 
occurring within the cell, chemical diffusion through the 
cell membrane, and intercellular electron transfer from the 
neighboring cell i — 1 to cell i, and then to the neighboring 
cell i + 1. The evolution of is also influenced by the 

external state s^^(f) experienced by the cell. Let = 

tnCH,TlATP,<l^EXT^^^EXT^ 

time t, and f+ be the time instant immediately following t. 
We define the following processes which affect the evolution 
of (see Fig. [^; 

• Electron donor diffusion through the membrane, 

joining the internal electron carrier pool with 
rate Acjr(sj*^(f); s^^(f)). The state becomes 

= (jtlCH ^t^IATP-: QeXT’^EXT^ ’’ 

• Intercellular electron transfer from the neighboring cell 
z — 1, joining the high energy external membrane with rate 

(ill owing to the high transfer 
rate approximation, the electron joining the low energy 
external membrane of cell z — 1 is immediately transferred 
to the high energy external membrane of cell z), resulting in 

) — y^CH tPatpi Hextt 9ext ^J ’ 

• Conventional ATP synthesis: this process involves the trans¬ 

fer of one electron from the internal electron carrier pool 
to the internal membrane with rate /icjr(f), m 
resulting in the synthesis of one unit of ATP. The electron 
then leaves the internal membrane and either follows the 
aerobic pathway (i.e., it is captured by an electron accep¬ 
tor), with rate /io( 7 T(s/^(f); s^^(f)), resulting in = 

{mcH — ^jPATP + ^tQ^extiH^i^xt)^ anaerobic one 

(i.e., it is collected in the low energy external membrane, also 
the high energy external membrane of cell z -I- 1), with rate 
XEXT{s['\t);s^^\t)), so that s^'^(t+) = {mcH-l,nATP + 
Ij Qext + ^ext)’ 

• Unconventional ATP synthesis: this process involves the 

transfer of one electron from the high energy exter¬ 
nal membrane to the internal membrane to synthesize 
one unit of ATP, with rate After¬ 

wards, the electron follows either the aerobic pathway 
with rate /zoc/T(s 7 *^(f); s^^(t)), resulting in = 

{mcH,nATP + ^,q'j^xT^9^EXT ~ 1)’ anaerobic 

one with rate (t)), so that = 

ijtlCHittATP + ^^^EXT ~ ^)’ 

• ATP consumption, with rate p.ATpi^^i'^ (^); resulting 

in Sj (F*") = ijncH ? tiATP 1) Qext^ Hext')' 

A. Information capacity of bacterial cables 

In 1^ , we have studied the information capacity of bacte¬ 
rial cables. By treating the bacterial cable as a communication 
medium, its capacity represents the maximum amount of 
information that can be transferred through the cable, and thus 
sets the limit on the rate at which cells can communicate via 
electron transfer. 

Importantly, as shown in | [M) and in the following analysis, 
communication over a bacterial cable entails two conflicting 
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Fig. 4. Communication system over a bacterial cable. "ENC"=encoder, 
’'DEC"=decoder (64| . 

factors; 1) achieving high instantaneous information rate; 2) 
inducing the bacterial cable to operate in information efficient 
states. In fact, the input electron signal affects the energetic 
state of the cells along the cable via the electron transport 
chain, and thus affects their ability to relay electrons and to 
transfer information. 

Indeed, our analysis reveals that the size of the ATP queues 
along the cable is extremely important to the health of the 
cable and to its capacity to transfer information (or electrons). 
Thus, with a greedy approach (denoted as myopic signaling 
and analyzed in Sec. |III-A3| l, which ignores these biological 
constraints, the overall achievable rate is measurably reduced 
compared to the optimal capacity achieving scheme, analyzed 


in Sec. III-A2 Thus, the importance and impact of including 
the unique constraints and dynamics of bacterial cables is 
underscored. 

1) Queuing model: The model m presumes that the 
channel state is given by the interconnection of the state 
of each cell in the cable, leading to high dimensionality. 
As shown in Fig. the bacterial cable contains N cells. 
Letting Sceii be the state space of the internal state of each 
cell, then the overall state space of the cable is 5 = 
which grows exponentially with the bacterial cable length 
N. Indeed, one of the advantages of our proposed queuing 
model is its amenability to complexity reduction. Herein, we 
propose an abstraction of m by treating the bacterial cable 
as a black box, which captures only the global effects on 
the electron transfer efficiency of the cable, resulting from 
the local interactions and cells’ dynamics presumed by ID- 
Specifically, we let E{t) be the number of electrons carried in 
the cable, i.e., the sum of the number of electrons carried in 
the external membrane of each cell, which participate in the 
electron transport chain to produce ATP for the cellj^The state 
space for this approximate model is £" = {0,1,..., i?niax}, 
where E’max is the electron carrying capacity of the bacterial 
cable. Letting L'max^ be the electron carrying capacity of a 
single cell, we have that i?max = N ■ Em&x\ which scales 
linearly with the cable length, rather than exponentially. 

The rationale behind this approximate model is as follows: 
when E{t) is large, the large number of electrons in the 
cable can sustain a large ATP production rate, so that the 


^Our model in |64) studies a more general model with leakage and 
interference. 


ATP pools of the cells are full and the cable is clogged. 
When this happens, the electron relaying capability of the 
bacterial cable is reduced, so that only a portion of the input 
electrons can be relayed. On the other hand, when E{t) is 
small, a weak electron flow occurs along the cable, so that the 
ATP pools are almost empty, the cells are energy-deprived, 
and the cable can thus sustain a large input electron flow to 
recharge the ATP pools. We thus define the clogging function 
a{E{f)) € [0,1], which represents the fraction of electrons 
which can be relayed at the input of the cable. This model 
is in contrast to classical relay channels over cascades, where 
the performance is dictated by the worst link rather than by a 
clogging function IB- 

The communication system includes an encoder, which 
maps the message m G {1,2,...,M} to a desired electron 
intensity A(f), taking binary values X{t) G {Amin, Amax}, f S 
[0, r]j^ where T is the codeword duration, and Amin > 0 and 
Amax > Amin are, respectively, the minimum and maximum 
electron intensities allowed into the cable. 

The channel state E{t) evolves in a stochastic fashion as 
a result of electrons randomly entering and exiting the cable. 
Electrons enter the bacterial cable following a Poisson process 
with rate Xin{t)=a{E{t))\{t), resulting in the state E{t) to 
increase by one unit, where a{i) G [0,1] is the clogging 
function and X{t) is the desired electron intensity, to be 
designed. Electrons exit the cable following a Poisson process 
with rate p{E{t)), resulting in the state E{t) to decrease 
by one unit. In general, a{i) is a non-increasing function of 
i, with a(i?max)=0, a(0)=l, a(i)>0,Vi<i?max- In fact, the 
larger the amount of electrons carried by the cable E{t), the 
more severe ATP saturation and electron clogging, and thus 
the smaller a{E{t)). 

2) Capacity analysis: Using results on the capacity of 
hnite-state Markov channels | |65| , we have proved that the 
capacity of bacterial cables is given by 

Era&x. 

1= max Tri{i)I{X{{);i), (4) 

^ [.^min Amax] ^_ q 

where 


• A(i) is the average desired input electron intensity in state 
E{f) = i, dehned as 

A(i)^E[A(f)|L;(f)=z]; (5) 

• I{X{i)]i) is the instantaneous mutual information rate 
in state E{t) = i with expected desired input electron 
intensity X{i), given by 

I{x;i) =a(i)a:log 2 

+ Aminlofof^); (6) 

^max '^min \ '^max / 

• is the asymptotic steady-state distribution of the 
bacterial cable, induced by the desired input electron 


^Herein, we restrict the input signal to take binary values, since this choice 
is optimal |64| . 
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Fig. 5. Achievable information rate as a function of the minimum clogging 
state value, amin? for the MP and OPT input distributions |64|. 


intensity A(-), given by 


= n 

k=0 


a{k)X{k) 
fi{k + 1 ) 




(7) 


and 7r)^(0) is obtained via normalization 7 '‘a(*) = !)■ 

The capacity optimization problem in Eq. 0 is a Markov 
decision process with state E{t), action A(i) in state 
E{t) = i, which generates the binary intensity signal, and 
reward function /(A(i);i) in state E{t) = i, and can thus 
be solved efficiently using standard optimization algorithms, 
e.g., policy iteration (see ||66)). The following trade-off arises; 
1) the optimal input signal should, on the one hand, achieve 
high instantaneous information rate, i.e., it should maximize 
I{X{E(t))-, E{t)) at each time instant t; 2) on the other 
hand, it should induce an "optimal" steady-state distribution 
of the cable, such that those states characterized by less 
severe clogging and where the transmission of information is 
more favorable are visited more frequently. These two goals 
are in tension. In fact, the instantaneous information rate is 
maximum in states with large clogging state a{i) ~ 1, i.e., 
when E{t) is small and the bacterial cable is deprived of 
electrons. Visits to these states are achieved more frequently by 
choosing A(t) = Amin with probability one. However, under a 
deterministic input distribution the information rate is zero. 

3) Myopic signaling: While the optimal signaling A*(z) 
balances this tension by giving up part of the instantaneous 
information rate in favor of better states in the future, the 
myopic signal greedily maximizes the instantaneous informa¬ 
tion rate, without considering its impact on the steady-state 
distribution of the cable. This is defined, Vz G £, as 


Amp(*) — argmax J(a;;z) = ^^ 

[Amin ; Amax] ^ 


, ( 8 ) 


and is constant with respect to the cable state E{t) = i, 
so that it does not require channel state information at the 
encoder. Importantly, the myopic signal neglects the steady- 
state behavior of the cable, and thus it tends to quickly 
recharge the ATP reserves of the cells, resulting in severe 
clogging of the cable. Thus, the myopic signaling may induce 
frequent visits to states characterized by small clogging state 


a{E{t)) ~ 0, i.e., when E{t) approaches i?max- 

Indeed, in we have shown that any input distribution 
X{i) larger than the myopic one XMp{i) is deleterious to 
the capacity for the following two reasons: 1) a lower in¬ 
stantaneous mutual information rate is achieved, compared 
to Amp(*) (by definition of the myopic distribution, which 
maximizes /(x; z)); 2) faster recharges of electrons within the 
cable are induced, resulting in frequent clogging of the cable, 
where the instantaneous information rate is small. 

4) Numerical results: We consider a cable with electron 
capacity E’max = 1000. The clogging state a{i) and output 
rate p(z) are given b}0 


a(z) = x(z < -Bmax) 


1 (1 CZmin) 


E„ 


(9) 


/i(z) = 0.6 -I- 0.8 


Er, 


In Fig. 1^ we note that the achievable information rate 
increases with amin under both the optimal signaling (OPT) 
and the myopic signaling (MP). This is because, as amin 
increases, both a(z) and the instantaneous information rate 
/(x;z) increase as well (see Eqs. (j^ and (|^). Intuitively, 
the larger a(z), the better the ability of the bacterial cable 
to transport electrons. 

OPT outperforms MP by ~ 9% for small values of amin- 
In fact, clogging is severe when the cable state approaches 
the maximum value £^max- Therefore, in order to achieve 
high instantaneous information rate, the state of the cable 
E{t) should be kept small. MP greedily maximizes the in¬ 
stantaneous information rate, but this action results in an 
unfavorable steady-state distribution, such that the cable is 
often in large queue states E{t) ~ E'maxj where clogging is 
severe and most electrons are dropped at the cable input. On 
the other hand, OPT gives up some instantaneous information 
rate in order to favor the occupancy of low queue states, 
where a(z) approaches one and the transfer of information is 
maximum. The performance degradation of MP with respect 
to OPT is less severe when amin—>1- In fact, in this case 
the instantaneous information rate /(x;z) is the same in all 
states (except T^max, where a(E;max)=0 and /(x; E;max)=0), 
hence optimizing the steady-state distribution of the cable is 
unimportant. Further discussion can be found in |[64). 


IV. Towards a Model of Quorum Sensing in a 

HOMOGENEOUS CELL COLONY 


In this paper, we propose stochastic queuing models as a 
general framework to capture signal interactions in microbial 
communities. We apply this general framework to two differ¬ 
ent microbial systems: electron transfer over bacterial cables 
and quorum sensing. In the previous section, we summarized 
a queuing model of electron transfer over bacterial cables 
and we have shown how complexity reduction can be 
achieved by a more compact queuing model, which represents 
the number of electrons carried in the cable by the state 
variable E{t) and defines a clogging function a{E{f)) over 


^The specific choices of a(i) and zi(i) have been discussed with Prof. 
M. Y. El-Naggar and S. Pirbadian, Department of Physics and Astronomy, 
University of Southern California, Los Angeles, USA. 
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I Virulence 
factors 


Fig. 6. Quorum sensing signaling system. 


its State space (64). Similarly, in this section, we propose a 
queuing model for the simplest quorum sensing signaling, for 
the case of a homogeneous cell colony which uses a single 
autoinducer-receptor pair, depicted in Figs. 1^ and 1^ This is 
the case, for instance, for Chromobacterium violaceum m- 
Similarly to the queuing model of electron transfer, this model 
can then be used as a baseline for further model reduction 
techniques. 

We consider a colony composed of N{t) identical cells at 
time t. We assume that these cells can duplicate, but do not 
die. In fact, only a small fraction of cells under favorable 


growth conditions are dead (< 0.01%) |68|, therefore cell 
death has a negligible influence on quorum sensing dynamics. 
We model cell duplication as a Poisson process with intensity 
p{n) per cell, function of the cell population size N{t) = n, 
i.e., each cell duplicates once every 1/p{n) hours, on average. 
Therefore, N{t) is a counting process with time-varying 
intensity N{t)p{N{t)). One model of p{n) is the logistic 
growth model p{n) = Pmax (1 - n/N^^^) ||^, j^, where 
Pmax is the maximum growth rate and iVmax is the maximum 
population size. This model presumes that growth slows down 
as the population increases, until N{t) = A%ax when growth 
stops. We denote the volume of a single cell as (j)ceU, so 
that the total cell volume of a colony of N{t) = n cells is 
Vceiiin) = n(f>ceU, and the total volume occupied by the colony 
(cell volume and extracellular environment) as Vtot(n). Here 
we consider both closed systems, with a hxed total volume 
Vtot(n) = constant, and open systems, in which the volume 
changes with the number of cells, Vtotin) = n(j)ex, where 
(j>ex is the volume occupied by each cell (cell volume and 
extracellular environment). 

The quorum sensing signaling system is characterized by 
different signals, represented in Fig. autoinducer molecules, 
produced by the synthases within each cell and released in the 
environment; receptors, located within each cell, which bind 
with the autoinducer molecules to form autoinducer-receptor 
complexes', these complexes then bind to the active DNA sites 
to produce further synthases, receptors and virulence factors. 
The dynamics of each of these signals and their interaction 
are characterized herein. 


A. Environment state 

Let A{t) be the number of autoinducer molecules in the 
system. Thus, the concentration of autoinducers in the volume 



Fig. 7. Queuing model of quorum sensing signaling system. 


occupied by the colony is 


VA{t) = 


m 

Vtot{N{t)y 


( 10 ) 


The dynamics of A{t) over time is affected by autoinducer 
production, degradation, receptor binding, and leakage. Leak¬ 
age accounts for the autoinducer molecules that leak out of the 
volume before being captured by the cells. We assume that the 
larger the colony, the smaller the leakage. In fact, the larger the 
colony of cells, the higher the likelihood that these molecules 
are detected before they diffuse out of the system. Addi¬ 
tionally, autoinducer molecules may chemically degrade thus 
becoming inactive. We thus let SA(n) be the rate of leakage 
and degradation for each autoinducer molecule as a function of 
the number of cells N{t)=n, where (5^(n)>5yi(n-|-1), Vn>l. 
An example is 


SA{n) = + 


^L,l 

1 +CL.2(n 


1 )’ 


( 11 ) 


where is the degradation rate of each autoinducer molecule, 
is the leakage rate when only one cell is present, and 2 
is the rate of leakage decay as the colony size increases. Given 
that there are A{t) autoinducer molecules in the system, the 
overall leakage rate is A(f)i5^(W(f)). 


Remark 1. Both 1 and 2 depend on the diffusion proper¬ 
ties of the medium and the spatial distribution of the cells GD- 
For a closed system without leakage, we have that 1 = 0, 
hence 5A{n) = ^D,Vn. However, (111 may not capture all 
possible scenarios. Other scenarios can be accommodated by 
an appropriate selection of SA{n) (possibly, different from 


([n]i). 

B. Cell state 


Each cell is described by three state variables, as represented 
in Fig. I?] For cell i, we let Ri (t) be the number of unbound 
receptors at time t, and Ci{t) be the number of autoinducer- 
receptor complexes within the cell. Each cell uses a protein 
called synthase to synthesize autoinducers. In turn, these au¬ 
toinducers are released outside the cell through its membrane. 
We let Si {t) be the number of synthases within cell i at time t. 
The state of cell i is thus given by (iij(f), Ci{t), Si{t)) at time 
t. No leakage of receptors, complexes or synthases outside of 
the cell occurs. 
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Each of the Ci (f) autoinducer-receptor complexes randomly 
binds to active sites in the DNA sequence, resulting in specific 
gene expression (see Fig. [^. We assume three active binding 
sites: 1 ) the first site contains the code to produce synthases, 
resulting in the increase of Si{t) by one unit; 2) the second 
site contains the code to produce receptors, resulting in the 
increase of Ri{t) by one unit; 3) the third site contains the 
code to produce the virulence factor. 

Remark 2. It is worth noting that, while in the case of 
Chromobacterium violaceum the receptor and the synthase 
are made from distinct DNA binding sites, in some cases 
both compounds are made from the same DNA binding 
site ilZD- This scenario can be incorporated by letting the 
first site contain the code to synthesize both synthases and 
receptors, resulting in the increase of both Si{t) and Ri{t) by 
one unit. Additionally, individual binding events may lead to 
the production of multiple synthases or receptors | |7^ . This 
scenario can be included by defining a random number of 
synthases or receptors produced at each binding event, with a 
given probability distribution. 

After the gene is expressed, the autoinducer-receptor com¬ 
plex unbinds from the DNA site. We assume that the binding¬ 
unbinding process is instantaneous, i.e., we neglect the amount 
of time that the complex occupies the DNA site. We model 
the process of autoinducer-receptor complexes binding to the 
DNA sites as Poisson processes, whose intensities depend on 
the affinity of the complex to a specific site. Thus, we let 
G {1; 2, 3} be the binding rate of each complex to site 
j (see Fig. [^. Additionally, gene expression occurs even in 
the absence of complexes, at basal rate eoj,j G {1,2,3}. 
Therefore, since there are Ci(t) autoinducer-receptor com¬ 
plexes within cell i, the gene located in site j is expressed 
with overall rate cqj + Ci{t)ec,j- 
C. State dynamics 

N{t) increases over time due to cell growth; A{t) increases 
over time as the synthases located within each cell produce 
autoinducers, and decreases over time as a result of autoinduc¬ 
ers binding to receptors, leakage and degradation (Sec. |IV-A] l; 
Ri{t) increases over time as the cell synthesizes receptors, 
and decreases over time as a result of autoinducers binding to 
receptors and receptor degradation; similarly, Ci{t) increases 
over time as the unbound receptors bind to the autoinducer 
signal in the environment, and decreases over time as a result 
of degradation and complex unbinding; finally. Si (t) increases 
over time as the cell produces synthases, and decreases over 
time as a result of degradation. We let 

\ CTOT(f) ^ C.it), (12) 


be the total amount of receptors, complexes and synthases of 
the cell colony, and 


VR{t) = 
vc{t) = 
Vsit) = 


-f^TOT(^) 

VcelliN{t))^ 

C'tot(0 

-S'tot(^) 

Vcell{N{t))^ 


( 13 ) 


be their concentrations, respectively]^ There are 
several processes that affect the dynamics of 
{N{t), A{t), RjoTit), Cjoj(t) ,Sjoj(t)), as detailed below. 

a) Cell duplication: Cell duplication occurs with rate 
p{N{t)) for each cell, so that the cell population N{t) 
increases by one unit with rate p{N{t))N{t). When one cell 
duplicates, its content is split randomly among the two cells. 
Therefore, if {Ri(t),Ci{t), Si{t)) is the state of cell i before 
its duplication, after the duplication (at time instant f+) 
there are two cells with state 

and respectively, where 

i?W(f+) + R<f\t+) = R^t), Cf^(f+) + = Qit), 

and + with probability distribution 

P((i?f ^(t+), Cf )(f+), )(t+))=(r, c, ait), 

^ ^ ^ ^ C^it) ^ ^ S,it) ^ ( 14 ) 

b) One autoinducer is created: Each synthase produces 
autoinducers with rate /3. Since there are 5TOT(i) synthases, 
the cell colony produces autoinducers with rate pSiaiit), 
resulting in the increase of A{t) by one unit. 

c) One autoinducer leaks or degrades: Each autoinducer 
leaks or degrades with rate 5A{N{t)). Therefore, autoinducer 
leakage and degradation occurs with rate (5^(A^(f))A(f), re¬ 
sulting in the decrease of A(f) by one unit. Degradation rates 
are dependent on auto-inducer type and environmental factors, 
as discussed in lig, 1 ^. 

d) One receptor is created: One receptor is created 
whenever one complex binds to the second active DNA site. 
This occurs with rate W(t)eo ,2 + C'TOT(i)ec ,2 across the whole 
cell colony, resulting in the increase of i?TOT(f) by one unit. 

e) One receptor degrades: Each receptor degrades with 
rate Sr. Therefore, receptor degradation occurs with rate 
^RRiaiit), resulting in the decrease of i?TOT(f) by one unit. 

f) One complex is created: Typically, binding of autoin¬ 
ducers to receptors to form autoinducer-receptor complexes 
does not occur until a threshold concentration of autoinducers 
is achieved IZ3- We let this threshold be pa,*- Thus, if 
r]A{i) < 'nA,±, then no binding of autoinducers to receptors 
occurs. When the threshold is exceeded, i.e., rjAit) > 'nA,±, 
then the binding of autoinducers to receptors is a Poisson 
process with intensity 7 [per unit of autoinducer and receptor 
concentrations, per unit volume, per hour]. Since the con¬ 
centrations of autoinducers and free receptors is 77^1 (t) and 
rjRit), respectively, and the total cellular volume where these 
reactions occur is Vceii{N{t)) = N{t)4>ceU, complexes form 
with rate 


Rtotit), -^( 0 ) 

- 'lVAit)VRit)x{VAit) > VA,th)VceuiN{t)) 
Ait)Rjcn:{t) 


= 7- 


Vtot{N{t)) 


-xiMt) > PA,±Vtot{N{t))), ( 15 ) 


^Note that receptors, complexes and synthases are located inside the cell 
rather then in the extracellular environment, therefore their concentrations are 
calculated with respect to the total cell volume Vceii{^{t))- 
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TABLE II 

Description of the queuing abstraction for the quorum sensing model. 


Queue 

Arrival rate 

Service rate 

Explanation 

N{t) 


0 

cell population size 

m 

/35TOT(t) -1- UcCxOTli) 

SA(N{t))A{t) + RTOT(i), N{t)) 

number of autoinducer molecules 

Rtot(() 

-l- (7roT(t)ec,2 -1- ucCtot(() 

^HRTOrlt) -1- Ac(A{t), i?TOT(t), Af(t)) 

number of receptors 

Crarit) 

Ac(A(t), RTorf), N{t)) 

ScCjctTit) -1- UcCxOTlt) 

number of complexes 

Sjarit) 

-1- C-Tarif^c,! 

SsSjorit) 

number of synthases 


so that A{t) and decrease by one unitj^and Ctot(0 

increases by one unit. 

g) One complex degrades: Each autoinducer-receptor 
complex degrades with rate 5c- Therefore, autoinducer- 
receptor complex degradation occurs with rate (5cC'TOT(f), 
resulting in the decrease of CTOT(f) by one unit. 

h) One complex unbinds: Each complex may randomly 
unbind, and the constituent autoinducer and receptor molecules 
become active again. This event occurs with rate vc for each 
complex. Therefore, unbinding of complexes occurs with rate 
vcCxm{t), resulting in the decrease of CTOT(i) by one unit 
and in the increase of A{t) and i?TOT(i) by one unit. 

i) One synthase is created: One synthase is created 
whenever one complex binds to the first active DNA site. 
This occurs with rate iV(f)eo,i +C'TOT(i)ec,i across the whole 
cell colony, resulting in the increase of S'TOT(i) by one unit. 
Since synthases create autoinducers and autoinducer-receptor 
complexes create new synthases, the production of autoinducer 
can be considered to be autocatalytic. 

j) One synthase degrades: Each synthase degrades with 
rate 6s- Therefore, synthase degradation occurs with rate 
SsSxcrrit), resulting in the decrease of S'TOT(f) by one unit. 

k) Virulence factor expression: Einally, the gene located 
in site 3 is expressed with intensity eg ,3 + Ci(t)ec ,3 in cell i, 
resulting in the overall expression with intensity -/V(t)eo ,3 + 
C'TOT(i)ec ,3 over the whole cell colony. 

D. State representation and state space reduction 

Let R{t) = , RN{t){t)\ be the vector of re¬ 

ceptor concentrations, C{t) = [C'i(f), ( 72 (f),..., C'Ar(t) (f)] be 
the vector of autoinducer-receptor complexes concentrations, 
and S{t) = [S'i(f), S' 2 (t),..., 5'Ar(t)(f)] be the vector of syn¬ 
thase concentrations. These vectors have time-varying length 
N{f), due to cell duplication. Then, the state of the system at 
time t is described by the tuple {N{t),A{t), R(t), C{t)^ S{t)). 
We notice that the state space grows exponentially with the 
cell population size, similarly to the model of electron transfer 
in Sec. However, the analysis of state dynamics above 
highlights that (7V(f), A(f), i?TOT(f), CtotW, *S'TOT(i)) is a 
Markov chain with dynamics described in Table |II] hence 
complexity reduction can be achieved. Indeed, this stochastic 
model and such Markov property can be exploited to infer 
the overall quorum sensing dynamics, rather than tracking the 
dynamics of each specific cell, i.e., 

Pt{n,a,r,c,s) ( 16 ) 

= P {N{t)=n, A{t)=a, i?TOT(i)=T, CjoT{t)=c, SjoTit)=s) - 


®In this paper, for simplicity, we assume that one receptor binds to one 
autoinducer to form one complex. However, in most systems, two receptors 
bind an autoinducer. 


V. Simulations Results and Experimental Data 


In this section, we present simulation results for the ho¬ 
mogeneous cell colony considered in Sec. [TV] The param¬ 
eters are listed in Table imi We consider two different se¬ 
tups: open system and closed system, both with initial state 
(V(0),A(0),i?TOT(0),CTOT(0),5TOT(0)) = (1,0,0,0,0) and 
the following parameters. 

a) Open system: Here, cells in an open system exist as a 
densely packed biofilm attached tho a surface. Autoinducers 
may leak, with parameters i = 5000 and ^ l ,2 = 0.1, and 
the cell colony grows in an open space. Cells grow tightly 
packed together, so that the extracellular environment takes 
only 10% of the cell volume = ^-^fceii)- The maximum 
population size is IVmax = oo, so that cell growth does not 
slow down. The concentration of autoinducers is computed as 


in (lOi, over the total volume Vtot(lV(f)) = N{t)(pex, whereas 
the concentration of receptors, complexes and synthases is 


computed as in (13i, over the cell volume Vceii{N{f)) = 
Nit)fceU- 

b) Closed system: Autoinducers do not leak i = 

2 = 0). The cell colony grows in a closed space of 
volume Vtot{n) = O.lnL. The maximum population size is 
-Vmax = 1000, so that, when growth is complete, the cell 
volume takes 1% of the total volume. The concentration of 
autoinducers is computed as in ( [T0| , over a constant total 
volume Vtot{N{t)) — O.lnL, whereas the concentration of 
receptors, complexes and synthases is computed as in ( [T3] l, 
over the cell volume Vceii{N{t)) = N{t)(j)ceih which grows 
with the population size. 

Open and closed systems occur both in nature and in 
experimental testbeds: microdroplets 0^® and well-plates 
are examples of closed systems, whereas cells on agar 
plates | |6^ , fTS) or in the presence of flow | |79| are examples 
of open systems. 

We use Gillespie’s exact stochastic simulation algorithm 
18^ to simulate the dynamics of the cell colony. In Eigs.|^a-e, 
we plot the time series of cell concentration, the concentration 
of autoinducers, receptors, complexes and synthases, respec¬ 
tively. Erom Eigs. |^d-e, we notice that the concentrations of 
autoinducer-receptor complexes and of synthases are small and 
approximately constant before 6 hours and 8 hours for the open 
and closed systems, respectively, and then grow quickly and 
steadily afterwards. We deduce that quorum sensing activa¬ 
tion occurs after approximately 6 hours in the open system 
and after 8 hours in the closed system. Indeed, as can be 
observed in Eig. [^b, the concentration of autoinducers grows 
slowly before quorum sensing activation, until reaching the 
critical concentration r]A,th = 21.4[nM]. Once this threshold is 
reached, complexes start to get formed with rate given by Eq. 
(151, thus explaining the steep increase observed in Eig. |^d. 















TABLE III 

Simulation parameters (nM denotes nanomolar concentration) 


II 


Parameter 

Explanation 

Value [all of them per hour] 

Reference 

pmax 

Maximum duplication rate 

1 [per cell] 


4^cell 

Cell volume 

1 [fL] 


/3 

Autoinducer generation rate 

18 [per synthase] 


75 



Autoinducer degradation rate 

0.01 [per autoinducer] 


m 


PA,th 

Autoinducer concentration threshold 

21.4 [nM] 

calculated from | /^ 


eo,2 

Basal receptor generation rate 

80 [per cell] 

calculated from 

75 


^C,2 

Activated receptor generation rate 

3 [per complex] 

calculated from 

E 



Receptor degradation rate 

12 [per receptor] 


81 
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Complex generation rate 

3.5 [per nM receptor & autoinducer concentrations, per fL] 


ST 


Sc 

Complex degradation rate 

1.4 [per complex] 


ST 


VC 

Unbinding of complexes 

60 [per complex] 


ST 


eo,i 

Basal synthase generation rate 

80 [per cell] 

calculated from 

75 


<SC,1 

Activated synthase generation rate 

3 [per complex] 

calculated from 

75 


Ss 

Synthase degradation rate 

1 [per synthase] 



In turn, the high concentration of complexes, which frequently 
bind to the DNA sites, contributes to the buildup of synthases 
in the system, thus explaining the steep increase observed in 
Fig. @e. Finally, each synthase produces autoinducers with 
rate /3, thus contributing to the buildup of autoinducers and 
sustaining the quorum sensing positive feedback loop. 

Quorum sensing activation is faster in the open system, 
since a much higher cell density is achieved (Fig. [^a), 
resulting in a steeper growth of autoinducer concentration, and 
thus earlier quorum sensing activation. In the open system the 
parameter that accounts for the loss of autoinducers to the 
surroundings will depend on the system, and further work is 
needed to determined if closed systems in general have slower 
activation kinetics. 


In Fig. 0a, we note that, when quorum sensing activation 
occurs, cell density is much smaller in the closed system than 
in the open system. This is because, in the closed system, 
autoinducers do not leak, and thus a lower concentration 
of cells is sufficient to reach the critical concentration of 
autoinducers for activation. In closed systems autoinducers are 
produced by dispersed cells and get diluted and retained in 
the entire volume of the system, whereas in open systems 
autoinducers are produced in concentrated regions of cells 
but escape in the surroundings. The local concentration of 
producers and autoinducer escape are key factors that dictate 
quorum sensing activation. Similar dependencies have been 
observed in other nonlinear systems |84|. 


In Fig. 0c, we note that, in the open system, the concentra¬ 
tion of receptors drops after quorum sensing activation. This 
is due to the high rate of autoinducers binding to receptors 
to form complexes, thus reducing the availability of free 
receptors. A similar behavior can be noticed in the closed 
system. However, the drop in this case is followed by a 
build up of receptors. In fact, the lower concentration of 
autoinducers in this case limits the production of complexes, 
resulting in excess free receptors. 

Finally, we provide experimental data for a closed system. 
The genes required for quorum sensing were placed into Es¬ 
cherichia coli using a plasmid constructed in | [85| . The genetic 
constructs place a green fluorescent protein under control of 
quorum sensing. Cells were grown in LB media at 37 °C. 
Growing cells are diluted such that cells remain in the early 
to mid exponential phase of growth, before quorum sensing is 


activated, for several generations before measurement. Since 
most proteins are stable for long periods of time, maintaining 
cultures at low density for several generations is needed to 
initialize the population of cells to a quorum sensing off state. 
Fluorescence was monitored using a 96 well-plate reader i fTT] 
and cell concentration was measured by absorbance of light 
at 600nm. 

As seen in Fig. 0e, the dynamics of the green fluores¬ 
cent protein regulated by quorum sensing in experiments is 
qualitatively similar to the predicted expression of synthase: 
they both exhibit the characteristic switch from a low to high 
production rate after quorum sensing is activated. However, 
green fluorescent protein production is immeasurable or zero 
prior to quorum sensing activation, which suggests that it has 
a lower basal level of gene expression than that assumed for 
synthase in our model. The basal level of expression can be 
independently tuned for each quorum sensing regulated gene. 
Genes that are part of the quorum sensing machinery, such as 
synthases and receptors, likely need a higher basal expression 
rate to prime quorum sensing activation. The consequences of 
differences in the basal and activated rates of expression for 
quorum sensing controlled genes warrants further investiga¬ 
tion. 

Fig. 0f plots the cell concentration over time and quorum 
sensing activation time. Note that the dynamics of cell concen¬ 
tration over time for the closed system (Fig. 0 a) are similar 
to experimental data (Fig. 0f), i.e., the cell population grows 
exponentially fast until reaching a certain maximum concen¬ 
tration, thus validating the logistic growth model employed. 
However, both the cell duplication rate and the maximum 
concentration exhibit different values in our model (pmax = 1 
[per cell per hour] and IVniax/Vtot(n) = 10^° [cells per 
mL], respectively) and in the experimental data (pmax = 0.65 
[per cell per hour] and Afmax/Vtot(ft) = 8 x 10® [cells per 
mL], respectively). Indeed, these parameters are variable and 
depend on specifics of each experiment, such as the cell type, 
composition of the media, and growth conditions. 

VI. Extensions 


Our framework can accommodate several aspects of quorum 
sensing observed in nature. 

A. Signal Integration 


Indeed, most quorum sensing signaling systems are more 
complex than the scenario considered in Sec. IV which 
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includes one receptor-autoinducer pair only. For instance, 
Pseudomonas aeruginosa contains multiple quorum sensing 
networks that work together to monitor changes in the lo¬ 
cal population of cells. Since the initial discovery of the 
LasI/LasR system, two other quorum sensing circuits in Pseu¬ 
domonas aeruginosa have been discovered and characterized, 
the Rhll/RhlR and PQS | | 8 ^ . Recently a fourth quorum 
sensing system has been identified, IQS p7) . Although each 
of these networks has a unique signal and receptor, these four 
systems are coupled in a feedback network in which each 
signaling pair regulates the production of both the synthase 
and receptor of at least one other network. Vibrio species 
have similar complexity in their quorum sensing networks 
188 ) , integrating information from multiple signals to direct 
regulatory decisions. 

B. Interference 

In addition to individual microbes producing multiple sig¬ 
nals, quorum sensing systems also mediate communication 
between different species. Many bacteria find themselves 
in diverse and crowded environments, and many of these 
neighboring species will also produce autoinducers. The type 
of autoinducer produced by many Vibrio and Pseudomonas 
species has been identified in more than 70 different species 
| [89) , and there are examples of signals or enzymes that inhibit 
quorum sensing activation | [90| . Within these multispecies 
communities, there is significant potential for interference with 
the process of quorum sensing activation, i.e., any process 
which alters the ability of cells to produce, exchange, rec¬ 
ognize, or respond to a quorum sensing autoinducer. Several 
mechanisms of interference that occur naturally have been 
identified. 

1) Signal synthesis inhibition: Although no examples of this 
interaction have been found in nature, a synthetic com¬ 
pound that prevented acyl-homoserine lactone (HSL) 
production by the Tofl synthase in the soil microbe 
Burkholderia glumae has been identified iD. This 
compound occupies the active site of the synthase, 
preventing autoinducer production. 

2) Destruction or chemical modification of the autoinducer: 

The common example of signal destruction is the AiiA 
enzyme isolated from Bacillus AiiA is an en¬ 

zyme that degrades a common type of autoinducer, 
thereby preventing the autoinducer from binding to the 
receptor. In another example, the oxidoreductase from 
Burkholderia is also capable of inactivating some autoin¬ 
ducers through chemical modification, thereby altering 
the specificity for receptor molecules. 

3) Competitive binding to the receptor: Some species have 
been shown to produce small molecules that interfere 
with autoinducer binding to the receptor, for example 
a furanone compound produced by the alga Delisea 
pulchra | [92) . In other cases the competition is from 
analogous autoinducers. There are many versions of 
the autoinducer HSL, with differing carbon tails, and 
each bacteria produces one or more types of HSL. 
For example. Pseudomonas aeruginosa produces the au¬ 
toinducers C4-HSL and 3oxo-C12-HSL, whereas Chro¬ 


mobacterium violaceum produces C 6 -HSL |93[ . The 
receptor of Chromobacterium violaceum has been found 
to be activated by C 6 -HSL, but deactivated or inhibited 
by HSL with longer carbon tails, such as 3oxo-C12-HSL 
| |94| . In this way, crosstalk between autoinducers pro¬ 
duced by neighboring bacteria can influence on signal 
transduction and quorum sensing activation. 

Our model in Sec. |IV]can be extended to include different 
types of interference. Consider an interfering signal with 
concentration rjjif) = I{t)/Vtot{N{t)), where I{t) is the total 
amount of interfering molecules in the volume of interest. The 
interfering signal is injected in the system with rate pi, and 
leaks with rate 5i{N{t)), which is a function of the colony 
size, similarly to autoinducers. We may consider the following 
models of interference: 

• Receptor inhibition: the interfering signal binds with the 

receptors, thus competing with the autoinducer signal; 
unlike the autoinducer-receptor complex, which induces 
specific gene expression by binding to the DNA sites, the 
complex formed by the interfering signal and the receptor 
becomes inactive; the binding of interfering molecules 
to receptors is a Poisson process with intensity 7 /^ 
[per unit of interfering signal and receptor concentra¬ 
tions, per unit volume, per hour]. Since the concentra¬ 
tions of the interfering signal and of free receptors is 
777 (t) and 7777 (f), respectively, the overall binding rate 
is rVi{' t)VR{t)^ceii{N{t)) , since these reactions occur 

inside the cells and the total cell volume is Vcea{N{f))-, 
correspondingly, both Rjcfrit) and I{t) decrease by one 
unit; 

• Synthase blocking: the interfering signal binds with 
the active site of the synthase, thus preventing autoin¬ 
ducer production; the binding of interfering molecules 
to synthases is a Poisson process with intensity 77 s 
[per unit of interfering signal and synthase concentra¬ 
tions, per unit volume, per hour]. Since the concen¬ 
trations of the interfering signal and of synthases is 
777 (t) and 77 s (f), respectively, the overall binding rate is 
liSili{t)vs{t)Vceii{N{t)); correspondingly, both S'TOT(f) 
and I{t) decrease by one unit; 

• Autoinducer degradation: the interfering signal destroys 
or chemically modifies the autoinducer signal, thus pre¬ 
venting it from binding to the receptor; the binding 
of interfering molecules to autoinducers is a Poisson 
process with intensity Sja [per unit of interfering signal 
and autoinducer concentrations, per unit volume, per 
hour]. Since the concentrations of the interfering signal 
and of autoinducers is 777 (f) and 77^1 (f), respectively, 
and these reactions occur over the volume Vtot{N(t)) 
occupied by the cell colony, the overall binding rate 
is SiAr]i{t)VA{t)Vtot{N{t))-, correspondingly, both A{t) 
and /(f) decrease by one unit; this type of interference 
can be also represented as an additional source of leakage 
for the autoinducers, so that their overall leakage rate is 
[6AiN{t)) -f SiAViit)]A{t). 

The model can be extended to include multiple interfering 
signals in a similar fashion. 
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VII. Conclusions 

In this position paper, we have explored how abstractions 
from communications, networking and information theory 
can play a role in understanding and modeling bacterial 
interactions. In particular, we have examined two forms of 
interactions in bacterial systems: electron transfer and quorum 
sensing. We have presented a stochastic queuing model of 
electron transfer for a single cell, and we have showed how the 
proposed single-cell model can be extended to bacterial cables, 
by allowing for electron transfer between neighboring cells. 
We have performed a capacity analysis for a bacterial cable, 
demonstrating that such model is amenable to complexity 
reduction. Moreover, we have provided a new queuing model 
for quorum sensing, which captures the dynamics of the 
quorum sensing signaling system, signal interactions, and cell 
duplication within a stochastic framework. We have shown that 
sufficient statistics can be identified, which allow a compact 
state space representation, and we have provided preliminary 
simulation results. Our investigation reveals that queuing mod¬ 
els effectively capture interactions in microbial communities 
and the inherent randomness exhibited by cell colonies in 
nature, and they are amenable to complexity reduction using 
methods based on statistical physics and wireless network 
design. Thus, they represent powerful predictive tools, which 
will aid the design of systems exploiting bacterial capabilities. 
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